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Abstract. We discuss the dimensional characterization of the solutions 
space of a formally integrable system of partial differential equations and 
provide certain formulas for calculations of these dimensional quantities. 

1. Introduction: what is the solutions space? 

Let £ be a system of partial differential equations (PDEs). We would like 
to discuss the dimensional characterization of its solutions space. 

However it is not agreed upon what should be called a solution. We can 
choose between global or local and even formal solutions or jet-solutions to 
a certain order. Hyperbolic systems hint us about shock waves as multiple- 
valued solutions and elliptic PDEs suggest generalized functions or sections. 

A choice of category, i.e. finitely differentiable C k , smooth C°° or analytic 
C w together with many others, plays a crucial role. For instance there are 
systems of PDEs that have solutions in one category, but lacks them in 
another (we can name the famous Lewy's example of a formally integrable 
PDE without smooth or analytic solutions, |L] ). 

In this paper we restrict to local or even formal solutions. The reason is 
lack of reasonable existence and uniqueness theorems (in the case of global 
solutions even for ODEs). In addition this helps to overcome difficulties 
with blow-ups and multi- values. 

If the category is analytic, then Cartan-Kahler theorem |Kaj guarantees 
local solutions of formally integrable equations |Goj and even predicts their 
quantity. We then measure it by certain dimension characteristics. 

If the category is smooth, formal integrability yields existence of solutions 
only if coupled with certain additional conditions (see for instance |Hoj ) . 
Thus it is easier in this case to turn to formal solutions, which in regular 
situations give the same dimension characteristics. With this vague idea let 
us call the space of solutions Sol(£). 

With this approach it is easy to impose a topology on the solutions space. 
However we shall encounter the situations, when the topological structure 
is non-uniform. 

To illustrate the above discussion, let's consider some model ODEs (in 
which case we possess existence and uniqueness theorem). The space of local 
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solutions for the ODE y' = y 2 is clearly one-dimensional, but the space of 
global solutions (continuous pieces until the blow-up) has two disconnected 
continuous pieces (solutions y = (a — x)^ 1 for a < or for a > 0) and a 
singular point (solution y = 0). Another example is the equation y' 2 +y 2 = 1, 
the solutions on(— e,e) of which form S , but the space of global solutions 
is M 1 (both united with two singular points in Sol(£)). 

We would like to observe the "biggest" piece of the space of Sol(£), so that 
in our dimensional count we ignore isolated and special solutions or their 
families and take those of connected components, that have more parameters 
in. 

It will be precisely the number of parameters, on which a general solution 
depends, that we count as a dimensional characteristic. Let us discuss the 
general idea how to count it and then give more specified definitions. 

Note that in this paper we consider only (over)determined systems of 
PDEs. Most results will work for underdetermined systems, but we are not 
concerned with them. 

2. Understanding dimension of the solutions space 

Let us treat at first the case of linear PDEs systems (the method can 
be transferred to non-linear case). We consider formal solutions and thus 
assume the system of PDEs £ is formally integrable. We also assume the 
system £ = £^ is of pure order k, which shall be generalized later. 

Thus for some vector bundle it : E(ir) — > M we identify £ as a subbundle 
£ k C J k {n) (see j^ lHollKLV^ and let £ t C J 1 (it) be its (I - fc)-th prolonga- 
tions, I > k. Then the fibres ££° C J^°(tt) at points x E M can be viewed as 
spaces of formal solutions of £ at x G M. To estimate size of we consider 
the spaces of linear functions on £\ jX , i.e. the space £* x . The projections 
nij-i : £i )X — > £i-i )X induce embeddings 7t*i_i ■ ^i_ lx £i x , and we have 
the projective limit 

Remark that £* is the module over all scalar valued differential operators 
on 7T, while the kernel of the natural projection J£°(7r)* — > £* can be viewed 
as the space of scalar valued differential operators on ir vanishing on the 
solutions of the PDEs system £ at the point x € M . Thus elements of £* 
are linear functions on the formal solutions E^. 

We would like to choose "coordinates" among them, which will estimate 
dimension of the formal solution space. To do this we consider the graded 
module associated with filtred module £*: 

g*(x)=@g!(x), 

where gi{x) are the symbols of the equation at x € M: 
gi(x) = £i, x /£i-i, x C S l T* <g> tt x 
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(we let £i = J (tt) for I < k), and reduce analysis of £ * to investigation of 
the symbolic module g*. 

This g* is the module over the symmetric algebra ST X M = ®S t (T x M) and 
its support Char^(£) C P C T*M is a complex projective variety consisting 
of complex characteristic vectors. The values tC p of the symbolic module g* 
at characteristic covectors p € C T* \ form a family of vector spaces over 
Char^(£), which we call characteristic sheaf. 

By the Noether normalization lemma ([E|) there is a subspace U C T X M 
such that the homogeneous coordinate ring ST X M/ Ann g* (x) of Chsx x (£) 
is a finitely generated module over SU. It follows that g*(x) is a finitely 
generated module over SU too. 

If g*(x) is a Cohen-Macaulay module (see [Ej, but we recall the definition 
later in a more general situation, then g*{x) is a free SJJ-module (we called 
the respective PDEs systems £ Cohen-Macaulay in IKL2I and discussed their 
corresponding reduction). 

Let a be the rank of this module, and p = dimC/. By the above discussion 
these numbers can be naturally called formal functional rank and formal 
functional dimension of the solutions space ££° at the point x € M, because 
they describe on how many functions of how many variables a general jet- 
solution formally depends (we shall omit the word "formally" later), or how 
many "coordinates" from £* should be fixed to get a formal solution. 

If the symbolic module is not Cohen-Macaulay, the module g*{x) over SU 
is not free, but finitely generated and supported on P C U*. Let ¥(U) be the 
field of homogeneous functions P/Q, where P,Q€ SU, Q / 0, considered 
as polynomials on U*. Thus ¥(U) is the field of meromorphic (rational) 
functions on U*. 

Consider ¥(U) (g) g*(x) as a vector space over ¥(U). Keeping the same 
definition for a, let us call the dimension of this vector space p formal rank 
of £ at the point x G M. 

It is clear that for Cohen-Macaulay systems the two notions coincide. 
However since g*(x) over SU is not free, we would like to give more numbers 
to characterize the symbolic module. 

Let us choose a base e 1; . . . , e r of ¥(U) ® g*(x) such that e%, . . . , e r are 
homogeneous elements of g*(x) and denote by Ti C g*(x) the SCZ-submodule 
generated by this base. It is easy to check that T\ is a free S'C/-module. For 
the quotient module M\ = g*(x)/T\ we have the following property: 

Ann/i ^ in SU, for any h £ M\. 

Therefore Ann Mi ^ and the support Hi of M\ is a proper projective 
variety in F C U*. 

We apply the Noether normalization lemma to Si, we get a subspace 
U\ C U , such that M\ is a finitely generated module over SU\. Its rank will 
be the next number p\ and we also get a\ = dim U\ , which we can call the 
next formal rank and formal dimension. 
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Applying this procedure several more times we get a sequence of varieties 
Ei and numbers (pi, which depends, in general, on the choice of the flag 
U D Ux D U2 D ■ ■ ■ and the submodules Tj of S77j_i. 

Thus we resolve our symbolic module via the exact 3-sequences 

-> Ti -» g* -> Mi -» over SZ7, 

-»• T 2 -»■ Mi -»• M 2 -► over SI/i, 

(with SuppMj = SuppTj+i ^ SuppMj+i) etc. 



3. Cartan numbers 

In Cartan's study of PDEs systems £ (basically viewed as exterior dif- 
ferential systems in this approach) he constructed a sequence of numbers 
Si, which are basic for his involutivity test. These numbers depend on the 
flag of subspaces one chooses for investigation of the system and so have no 
invariant meaning. 

The classical formulation is that a general solution depends on s p functions 
of p variables, s p _i functions of (p — 1) variables, . . . , si func tions of 1 
variable and sq constants (we adopt here the notations from |BCG 3 ; in 



Cartan's notations O 

we should rather write 5^, Sp -\- Sp—\, Sp -\- s^— 1 -\- Sp—2 
etc). However as Cartan notices just after the formulation [Uj, this statement 
has only a calculational meaning. 

Nevertheless two numbers are absolute invariants and play an important 
role. These are Cartan genre, i.e. the maximal number p such s p 7^ 0, but 
s p+ i = 0, and Cartan integer a = s p . 

As a result of Cartan's test a general solution depends on a functions of 
p variables (and some number of functions of lower number of variables, but 
this number can vary depending on a way we parametrize the solutions). 
Here general solution is a local analytic solution obtained as a result of 
application of Cartan-Kahler (or Cauchy-Kovalevskaya) theorem and thus 
being parametrized by the Cauchy data. 

Hence we can think of p as of functional dimension and of a as of func- 
tional rank of the solutions space Sol(f ). In fact, we adopt this terminology 
further on in the paper, because as was shown in the previous section it 
correctly reflects the situation. 

These numbers can be computed via the characteristic variety. If the 
characteristic sheaf over Char c (£) has fibers of dimension k, then 

p = dimChar c (£) + 1, a = k ■ degChar c (£). 

The first formula is a part of Hilbert-Serre theorem (0), while the second 
is more complicated. Actually Cartan integer a was calculated in 
in general situation and the formula is as follows. 

Let Char c (g) = U e E e be the decomposition of the characteristic variety 
into irreducible components and d e = dim/Cj; for a generic point x E S e . 
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Then 

a = ^2<l e • deg£ e . 

It is important that these numbers coincide with the functional dimension 
and rank of the previous section. Moreover the sequence of Cartan numbers 
Si is related to the sequence (pi,<Ji) of the previous section. 

This can be seen from the general approach of the next and following 
sections, which treat the case of systems £ of PDEs of different orders (we 
though make presentation for the symbolic systems, with interpretation for 
general systems being well-known [^} IKLV| |KL2l ). 

4. Symbolic systems 

Consider a vector space T of dimension n (tangent space to the set of inde- 
pendent variables, substitute to T X M) and a vector space ./V of dimension m 
(tangent space to the set of dependent variables, substitute to tt x = tt~ 1 (x)). 

Spencer 5-complex is de Rham complex of polynomial iV-valued differen- 
tial forms on T: 

-> S k T* ® N ^ S^T* <g> N ® T* 4 • • • 4 s k - n T* ® N <g> A n T* 0, 

where S l T* = for i < 0. Denote by 

6 V = i v o 8 : S k+1 T* ® N S k T* ® N 

the differentiation along the vector v € T. 

The Z-th prolongation of a subspace h C S k T* ® iV is 

={p£ S k+l T*®N : fi^ ... tf Wl p e fcVui, ... = ^T*®/inS fc+ 'T*®iV. 

Definition. ^4 sequence of subspaces gt- C S^T* ® iV, > 0, wrai/i go = N 
and gk C is called a symbolic system. 

If a system of PDEs £ is given as F\ = 0, . . . , F r = 0, where Fi are scalar 
PDEs on M, then T = TM, N ~ M m and the system g C ST* ® iV is 
given as /i = 0, . . . , / r = 0, where = cr(Fj) are symbols of the differential 
operators at the considered point (or jet for non-linear PDEs). 

With every such a system we associate its Spencer 5-complex of order k: 

0^g k ^ 9k-l <8> T* -^-> g fc _ 2 ® A 2 T* - . . . 

- 9i <g> A^T* -A, • • • -!+ 5fc _ n g, A n T* -> 0. 

Definition. The cohomology group at the term gi <8> A^T* is denoted by 
H l ' J (g) and is called the Spencer 5-cohomology of g. 

Note that g k = S k T* ® N for < k < r and the first number r = r m ; n (g), 
where the equality is violated is called the minimal order of the system. 
Actually the system has several orders: 

ord(g) = {k £ Z + \ g k ^ gfij. 
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Multiplicity of an order r is: 

m(r) = dimgj^/g,. = dhnH r ~ 1 ' 1 (g). 

Hilbert basis theorem implies finiteness of the set of orders (counted with 
multiplicities) : 

codim(g) := dim if*' 1 (g) = \^m(r) < oo. 

Starting from the maximal order of the system k = r max we have: g^+i = 
(l) 

al- 
ii we dualize the above construction over R, then Spencer ^-differential 

transforms to a homomorphism over the algebra of polynomials ST and 

9* = ®i9i becomes an ST-module. This module is called a symbolic module 

and it plays an important role in understanding PDEs. 

In particular, characteristic variety Char (g) C P C T* is defined as the 

support of this module Supp(g*) = {\p\ : (g*) p / 0} and the characteristic 

sheaf fC over it is the family of vector spaces, which at the point p G Char (g) 

equals the value of the module at this point JC P = g* jp ■ g* . For more 

geometric description see {HI lKLV| IKL2I . 

5. Commutative algebra approach 

We will study only local solutions of a system of PDEs £, which we 
consider in such a neighborhood that type of the symbolic system does not 
change from point to point (on equation) in the sense that dimensions of gk, 
of the characteristic variety Char (g) and of the fibers of IC are the same. 

It should be noted that if a system £ is not formally integrable and £' is 
obtained from it by the prolongation-projection method [KllM2| IKL2I , then 
the numbers p, a change in this process, i.e. either the functional dimension 
or the functional rank decrease. Thus from now on we suppose the system 
£ is formally integrable. 

The numbers p, a can be described using the methods of commutative 
algebra. Recall (HM1) that by Hilbert-S erre theorem the sum 

i<k 

behaves as a polynomial in k for sufficiently large k. This polynomial is 
called the Hilbert polynomial of the symbolic module g* corresponding to 
£ and we denote it by Pe(z). If p = deg Pe(z) and a = P^ p \z), then the 
highest term of this polynomial is 

P e (z) = az p + ... 

(see [H] for the related statements in algebraic geometry, the interpretation 
for PDEs is straightforward). 

A powerful method to calculate the Hilbert polynomial is resolution of 
a module. In our case a resolution of the symbolic module g* exists and 
it can be expressed via the Spencer <5-cohomology. Indeed, the Spencer 
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cohomology of the symbolic system g is IR-dual to the Koszul homology of 
the module g* and for algebraic situation this resolution was found in |Gr| . 
It has the form: 

® q H q - n ' n (g) ® ^ (BqH^+^-^g) ® ^ ... 

-> ©gF 9 - 1 ' 1 ^) ^ ® q H q '°(g) ® sHd ^ 5 * -> 0, 

where S^ -9 ' is the polynomial algebra on T^M with grading shifted by q, 
i.e. 5] 9 ' = S l ~ q T x M, and the maps have degree 0. 

Thus denoting W = dim H^(g) and r Q = dim S a TM = ( Q+ ™ _1 ) we 
have: 

dim 5i = Y,{h qfl n- q ~ h^Ti-g-i + h"' 2 Ti- q - 2 -■■■ + (-l) n h«> n Ti- q - n ). 

1 

Let also jp = Yla<f3 T a = dim JyM = (^ n ) be the dimension of the fiber 

of the vertical jets JyM, i.e. the fiber of the jet space J@M over M. Thus 
we calculate 

^dim 9i = J2i hq '°jk- g ~ h q > l j k - g ^ + h q > 2 j k _ q _ 2 - • • • ± h q > n j k _ q _ n ). 
i<k q 

Finally we deduce the formula for Hilbert polynomial of the symbolic 
module g* 



P £ (z) = ^(> o rf n ) - h"' 1 ('-^- 1 ) + 
q 

+ h q > 2 { z - q+ n n - 2 ) + (-l) n h q > n ( z - q )). 



Here 

n k ) = ^{z + l).{z + 2).--(z + k). 

Denote Sj(ki, . . . , k n ) = Yl hi"'ky the j-th symmetric polynomial 

ii<— <ij 

and let also 

n! 

Thus 

n + 1 n _ (n+l)(3n + 2) „ n _ n(n + l) 2 
~ ' 1 ~ 2 ~ 4^3! ' 3 ~ 2-4! ' 

_ (n+l)(15n 3 + 15n 2 -10n-8) 

^ = 48^1 etC ' 

If we decompose 
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then we get the expression for the Hilbert polynomial 



i,j,q V Jl ' k=0 



where 



k-j 



3=0 «.« y Ji- 

tS. Calculations for the Solutions space 

We are going to compute the dimensional characteristics of two important 
classes of PDEs. 

Involutive systems. These are such symbolic syste ms g = {gk} that all 
subspaces are involutive in the sense of Cartan (UJ lBCG 3 j (this defini- 
tion for the symbolic systems of different orders was introduced in |KLs| ). 
Thanks to Serre's contribution |GS| we can reformulate this via Spencer 
cohomology as follows. 

Denote by g^ the symbolic system generated by all differential corollaries 
of the system deduced from the order k: 



I*) 



S l T* <g> N, for % < k; 



9i "1 g^ k \ ioii>k. 

Then the system g is involutive iff H i,3 \g^) = for alH > k (this condition 
has to be checked for k £ ord(<?) only), see |KLs| . 

In particular, H l ' J (g) = for i £ ord(<?) — 1, ^ (0,0), and the 

resolution for the symbolic module g* as well as the formula for the Hilbert 
polynomial of 8 become easier. 

Let us restrict for simplicity to the case of systems of PDEs 8 of pure first 
order. Then 

P £ (z) = " h^CllX 1 ) + ^ 2 { z+ n lf) -■■ 

= b x - — - + 627 -r. H V b . 

(n — lj! [n — 2y. 

Vanishing of the first coefficient 60 = is equivalent to vanishing of Euler 
characteristic for the Spencer 5-complex, x = Yli(~ l) 1 ^ = 0, and this is 
equivalent to the claim that not all the covectors from C T* \ are charac- 
teristic for the system g. 

The other numbers bi are given by the general formulas from the previous 
section, but they simplify in our case. For instance 

b 1 = n±k b - ^i-iyhPH = ^{-iy +1 i ■ h°'\ 

If codim Char c (f ) = n — p > 1, then b\ = and in fact then 6^ = for 
i < n — p, but b n _ p = a. 
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Theorem. 7/codimChar {£) = n—p, then the functional rank of the system 
equals 

^ (n-p)\ 
Proof. Indeed one successively calculate the coefficients using the formula 

i o=0 

and notes that equals to the displayed expression plus a linear combina- 
tion of bk-i, . . . , 6q. The claim follows. □ 

One can extend the above formula for general involutive system and thus 
compute the functional dimension and functional rank of the solutions space 
(some interesting calculations can be found in classical works |Jl 

Cohen-Macaulay systems. A symbolic system g (and the respective 
PDEs system £ ) is called Cohen-Macaulay ( IKL2I ) if the corresponding sym- 
bolic module g* is Cohen-Macaulay, i.e. (see |M11 IE| for details) 

dim<7* = depth g* . 

Consider an important partial case (we formulate the definition only for 
symbolic systems; PDEs are treated in IKL4I ): 

Definition. A symbolic system g C ST* <g> N (n = dimT, m = dimN) of 
codim(<7) = r is called a generalized complete intersection if 

• m < r < n + m; 

• codimc Char (g) = r — m + 1; 

• dim/C x = 1 Vx G Char c (#) C P C T* . 

Formal integrability of such systems are given by the compatibility con- 
ditions expressed via brackets (for scalar systems |KLj| IKL3I ) or multi- 
brackets (for vector systems IKL^p . In this case we can calculate Cartan 
genre and integer directly. 

Theorem. Let £ be a system of generalized complete intersection type and 
suppose it is formally integrable. Then the functional dimension of Sol(£) 
is 

p = m + n — r — 1 

and the functional rank is 

o~ = SV— m+i {ki , • • • , k r ) = ^ k{ 1 • ■ ■ ki r _ m+1 , 

i\<—<i r - m -\-l 

the l-th symmetric polynomial of the orders k\, . . . , k r of the system. 

Note that if the last requirement in the definition of generalized com- 
plete intersection is changed to dim IC X = d everywhere on the characteristic 
variety, then the functional rank will be multiplied by d: 

o~ = d ■ S r — m -\-i(ki, . . . , k r ). 
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However the formal integrability criterion for generalized complete intersec- 
tions is proved in IKL4I under assumption that d = 1. 

Proof. We shall consider the case of a system g of a pure order: k\ = ■ ■ ■ = 
k r = k, ki £ ord(<7). The case of different orders is similar and will appear 
elsewhere. 

The formula for functional dimension p follows directly from the definition 
of generalized complete intersection. Let's calculate a. 

We can use interpretation of the Cartan integer a from SJH1 Recall that 
characteristic variety Char c (g) is the locus of the characteristic ideal 1(g) = 
Ann(g), which the the annihilator of g* in ST. 

Since the module is represented by the matrix with polynomial entries 
(each differential operator Aj giving a PDEs system £ is a column Ay, 
1 < i < r, 1 < j < m; so that their union is a m x r matrix M(A)), 
its annihilator is given by the zero Fitting ideal (in fact, here we use the 
condition on grade of the ideal: depth Ann(g) = r — m + 1, which follows 
from the conditions of the above definition). 

This ideal Fitto(<?) is generated by all determinants of m x m minors 
of the corresponding to M(A) matrix of symbols M(cja). These minors 
are determined by a choice of m from r columns, so that there are Q") 
determinants and each is a polynomial of degree k r ~ m+1 . 

However not all the minors are required to determine Char c (g) and this is 
manifested by the fact, that we sum (-JLi) degrees k r ~ m+1 to get the func- 
tional rank a. The easiest way to explain this is via the Hilbert polynomial 
of the symbolic module g* . 

This can be calculated since under the assumption of generalized complete 
intersection g* possesses a resolution in the form of Buchsbaum-Rim complex 
(see |KL71 ): 

-> S^-m-ly* ^ K rjj ^ S r-m-2y* ^ ^-ljj ^ 

► A m+1 C7 -» U -» F -» 5 * -» 0, 

where V ~ 5T ® iV* (recall that dimiV = m and 5 C ST* (g> N) and 
C7 = ST ■ ■ ■ © ST. Star ★ means dualization over ST and the tensor 

V v ' 

r terms 
products are over ST as well. 

Now the claim follows from the detailed investigation of degrees of the 

homomorphisms in the above exact sequence. To see this we suppose at 

first that r = m + n — 1 and use the following assertion. 

Lemma. The following combinatorial formula holds: 
/n+fe(n+m-l)\ _ ( n m _ ^m+fe(n+m-2)\ 

(n+m— 1\ i n 



m 

n-1 
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We would like to comment and interpret the sum on the left hand side of 
this formula. In our case the system is of finite type (g* has finite dimension 
as a vector space) and a = ^ dim gi (the sum is finite). 

Stabilization of the symbol occurs at the order i = ^ ki — 1 = k(n + m — 
1) — 1: gi = 0. So we prolong £ to the jets of order k(n + m — I) and the 
first term is just dim J^ n+m (M, N). 

The next term is due to the fact that £ C J k (M, N) is proper. It is given 
by r = n + m - 1 equations of order k, we which we differentiate up to 
k{n + m — 2) times along all coordinate directions (prolongation). 

There are relations between these derivatives. These are compatibility 
conditions (1-syzygy of g*), which appear in the form of multi-brackets 
IKL4I , in our case this bracket uses (m + l)-tuples of Aj. 

There are in turn relations among relations (2-syzygy of g*), which are 
identities between multi-brackets (these we call generalized Pliicker iden- 
tities, to appear soon), in our case these latter use (m + 2)-tuples of the 
defining operators Aj etc. 

Due to exact form of the relations (higher syzygies) we get factors ("'^^Y 2 ) 
in the summations formula of the lemma. 

In the case r < n + m — 1 we should perform a reduction, which is 
possible by Theorem A IKL2I . Then the functional dimension p grows, but 
the functional rank remains the same and the previous calculation works. □ 

7. Examples 

Here we show some examples demonstrating the above results. 

1. Intermediate integral of a system £ C J ir is such a system £ C J k Tt 
that k < k and £ C £( k ~ k ) (where £W is the i-th prolongation of the 
system). Since every solution to the system £ is a solution to £~( k ~ k ) we 
conclude: Whenever the functional dimension p > 0, we have p = p and 
a = a. 

Indeed the solutions of £~( k ~ k ) f orm a finite-dimensional parametric family, 
such that solutions of £ appear for some fixed values of parameters (because 
we differentiate with respect to all variables to obtain the prolongation). 
Thus the number of functions of p > variables, on which a general solution 
depends, will not be altered. 

2. If the PDEs system £ is underdetermined, then p = n and cr > 1. 
Indeed, a is precisely the under-determinacy degree, i.e. the minimal number 
of unknown functions that should be arbitrarily fixed to get a determined 
system. We assume we can do it to get a formally integrable system. If 
underdetermined system is not formally integrable, compatibility conditions 
can turn it into determined or over-determined and then decrease p and 
change a. 

A nice illustration is the Hilbert-Cartan system 



z\x) = (/(*)) 2 . 
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It has p = 1, a = 1. But even though a general solution depends on one 
function of one variable, it cannot be represented in terms of a function and 
its derivatives only (Hilbert's theorem). 

3. As we noticed earlier the similar situation happens to overdeter- 
mined system: If £ is not formally integrable, and £ is obtained from £ 
by prolongation-projection technique (sometimes it is said that £ is the in- 
volutive form of £, but this is not true, only a certain prolongation of £ is), 
then p < p or \p = p and a < cr]. Indeed, supplement of additional equations 
shrinks the solution space. 

For instance if we consider two second-order scalar differential equations 
on the plane 

F(x, y, u(x, y), Du(x, y), D 2 u(x, y)) = 0, 
G(x, y, u(x, y), Du(x, y), D 2 u(x, y)) = 0, 

such that F and G have no common complex characteristics, then the com- 
patibility condition of this system £ can be expressed via the Mayer bracket 
( IKLTI ): H = [F,G] £ . UH = 0, then p = 0, a = A. If H + 0, then p = and 
a < 3, the equality being given by the Frobenius condition for the system 
£ = {F = 0,G = 0,H = 0}. 

If the system has one common characteristic and is compatible, we have: 
p = 1, a = 1. Pairs of such systems are basic examples of Darboux integra- 
bility. 

4. Evolutionary equations ut = L[u] provide interesting examples, which 
usually "contradict" the theory. Consider for instance the heat equation 

1H — u X x- 

It is formally integrable and analytic. We can try to specify the initial 
condition u\t=o = an d then solve the Cauchy problem, so that we get 
p = 1, a = 1. On the other hand we can let u\ x= q = *po(t), u x \ x= o = ipi(t) 
and then get p = 1, a = 2. 

If we calculate the numbers using our definitions of functional dimension 
and functional rank (for instance, via Hilbert polynomial), it turns out that 
the second approach is correct. Indeed with the first idea we come into 
trouble with certain Cauchy data: Let, for instance, (p(x) = (1 — x) , 
which is an analytic function around the origin. Then the analytic solution 
should have the series 

1 2 t 4! t 2 | (2n)! t n 

U{t,X ' ^ T^ + I(l-x) 3 + 2!(l-x) 5 + "' + n\ (l-x) 2 ™+ 1+ "' 

which diverges everywhere outside t = 0. The reason why the second ap- 
proach provides no problem is because the line {x = 0} is non-characteristic 
and we can solve our first order PDE by the classical method of Cauchy 
characteristics. 
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Remark however that in the standard courses of mathematical physics the 
heat equation is solved with the first approach (by Fourier method). How 
is it possible? 

Explanation is that we solve the heat equation then only for positive time 
t > 0. Doing the same method in negative direction blows up the solutions 
immediately (heat goes rapidly to equilibrium, but we cannot predict even 
closest past)! We here are interested in the solutions, which exist in an open 
neighborhood of the origin (like in Cauchy-Kovalevskaya theorem), and this 
contradicts the first approach. 

5. Similar problems arise with Cauchy problems in other PDEs systems: 
one usually applies reduction or fixes gauge, but this can change dimensional 
characteristics. 

For instance, consider the Cauchy problem for the Einstein vacuum equa- 
tions, which is a system of 10 PDEs of 10 unknown functions. The system is 
over-underdetermined (i.e. it has compatibility conditions). In wave gauge 
|CB| its solution depends on several functions on a 3-dimensional space, 
which are subject to constraint equations, so that p = 2. On the other 
hand, the original Einstein system is invariants under diffeomorphisms and 
this yields p = 4. 

One should also be careful with Cauchy data in higher order, since then 
the definition of characteristics becomes more subtle, see |KLs| . 

6. Consider a system £ , which describes automorphisms of a given geo- 
metric structure. The corresponding symbolic system is g C ST* eg) T. The 
automorphism group has maximal dimension iff the system is formally inte- 
grable. Consider the examples, when the geometric structure is symplectic, 
complex or Riemannian (all these structures are of the first order). 

Let at first g be generated by g\ = sp(n) C T* (g> T. Our tangent space 
T = T X M is equipped with a symplectic structure uj, and we can identify 
T* ~ T and we get g\ = S 2 T* C T* <S> T*. The prolongations are g% = 
S i+1 T* C S' l T* ® T. 

The system is easily checked to be involutive and the only non-vanishing 
Spencer 5-cohomology groups are 



Then one checks that the Euler characteristic is x = 1 7^ an d so &o 0- 
Thus the functional dimension is p = n. Indeed the characteristic variety is 
P C T* because each non-zero covector p is characteristic: p 2 E g\ — S 2 T*. 
Next by a theorem from £jf)lone calculates the functional rank 



This result is easy to verify: an infinitesimal symplectic transformation has 
a generating function (Hamiltonian) and so it is determined by one function 
of n variables. 



H°*(g) = A i+1 T*. 
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If we turn to (almost) complex structures J on M, then g\ = gl(^,C) = 
T* <g)c T (space of C-linear endomorphisms of T) and the prolongations are 

9i = S l c T* ®c T. 

The characteristic variety is proper and one calculates that p = \, a = n. 
The system is again involutive. The second Spencer cohomology is 

H°' 2 (g)=A 2 c T* ® c T, 

which is the space of C-antilinear skew-symmetric (2, 1) tensors (Nijenhuis 
tensors) . 

The last example is the algebra of Riemannian isometries (i.e. T is 
equipped with a Riemannian structure) of a Riemannian metric q on M. 
The symbol is g\ = o(n) and the prolongations are zero gi = g$ = ■ ■ ■ = 0. 

This system is not involutive. For instance, 

H 1,2 (g) = Ker(S 2 A 2 T* -> A 4 T*) 

(the space of Riemannian curvatures) is non-zero (for n = dimT > 1). Since 
the system is of finite type, the characteristic variety is empty and p = 0. 
The general solution (isometry) depends on a = ( ra + 1 ) ra constants. 

We recall, that the above dimensional conclusions are correct if the system 
£ is integrable, otherwise the space Sol(£) shrinks. In the above examples 
this means: the form u is closed (with just non-degeneracy we have almost- 
symplectic manifold); the structure J is integrable (Nijenhuis tensor Nj 
vanishes); the manifold (M,q) has constant sectional curvature (so it is a 
spacial form). 
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